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Представлены быстрые алгоритмы морфологических преобразований бинарных изображений. Отличительной особенностью 
представленных алгоритмов является исключение операций поиска минимального или максимального элемента внутри струк- 
турирующего элемента, что существенно повышает скорость реализации данных алгоритмов. 


Использование методов морфологического 
анализа изображения для коррекции изображений, 
искаженных случайной некоррелированной поме- 
хой является одним из перспективных направле- 
ний в области анализа изображений [1—8]. Для 
оперативных систем обработки и анализа изобра- 
жений важную роль играют их скоростные характе- 
ристики. Поэтому построение эффективных, бы- 
стрых алгоритмов морфологических преобразова- 
ний изображения является актуальной задачей. 

Пусть К и 2 представляют соответственно мно- 
жества действительных и целых чисел, и пусть Е 
представляет ^/-мерное непрерывное пространство 
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К\ (с1=1,2,3,...,К) или дискретное пространство 2 Л . 
Тогда ^/-мерный сигнал можно представить как 
функцию области К Л (непрерывной) или 2 Л (дис- 
кретной), диапазон которой составляет либо К - 
при непрерывном изменении амплитуды, либо 2— 
при квантованном изменении амплитуды. Двоич- 
ные сигналы могут быть представлены с помощью 
множеств. Например, изображение в левой верх- 
ней части рис. 1 представляет двоичный сигнал, в 
котором область белого фона выражена через 0, а 
заштрихованная зона через 1. Ясно, что сигнал мо- 
жет быть также представлен множеством X точек, 
соответствующих заштрихованной зоне. 


Технические науки 


Двоичные изображения часто получают в ре- 
зультате пороговой селекции полутоновых изобра- 
жений, и пороговая селекция часто используется 
для того, чтобы представить полутоновые изобра- 
жения через двоичные сигналы, то есть с помощью 
множеств. Серра [1] использует представление дей- 
ствительной (/-мерной функции /*), где * означает 
^/-мерный вектор, с помощью ансамбля его (/-мер- 
ных пороговых множеств, определяемых как 
Т а(Л = {х-/(х)>а,-са<а<со}, 

где амплитуда а полностью перекрывает К или 2, в 
зависимости от того, имеет ли сигнал /непрерыв- 
ный или квантованный диапазон. Пороговые мно- 
жества линейно упорядочены, поскольку «</>=>- 
Г„(/)з /,(/), и позволяют однозначно восстановить 
сигнал, так как 

/(х) = та х{а:хе/(/)}, Ѵх. (1) 

Сигналы//*) несут ту же информацию, что и по- 
роговые множества Т а (х). Следовательно,/*) может 
быть восстановлено по ансамблю /(*), поскольку 

/(*) = шах {я : / а =1}, Ѵх. 



Исходный сигнал 


Эрозия Наращение 


Размыкание Замыкание 

Рис. 1 . Эрозия, наращение, размыкание и замыкание мно- 
жества точек X с помощью структурирующего элемен- 
та В. Штриховые области соответствуют внутренности 
множеств, сплошная черная линия обозначает грани- 
цы преобразуемых множеств, а пунктирная линия ~ 
границу исходного множества точек X 

Фундаментальными морфологическими опера- 
торами для множеств являются наращение 0 и 
эрозия О с помощью В, В — структурирующий эл- 
емент, представляющий собой компактное множе- 


ство малого размера и простой формы. Операция 
наращения всех пороговых множеств функции / с 
помощью одного и того же компактного множества 
дает множества, Т а (/)®В, которые являются поро- 
говыми множествами новой функции ]®В, назы- 
ваемой наращением функции / с помощью В. Эта 
новая функция может быть вычислена либо из ( 1 ) 
как (/05)(*)=тах{а:*б Т а (/)®В}, или из прямой эк- 
вивалентной формулы 

(/ © В)(х) = тах{/(* - у)} . 

уеВ 

Подобно этому операция эрозии всех порого- 
вых множеств функции / с помощью одного и того 
же множества В и суперпозиции всех выходных 
множеств посредством (1) дает новую функцию, 
называемую эрозией функции/с помощью множе- 
ства 5, которая может быть вычислена по эквива- 
лентной формуле 

(/ОВ)(х) = тіп{/(х + у)}. 

уеВ 

Другие операторы могут быть определены как 
комбинации эрозии и наращения. Например, два 
дополнительных фундаментальных оператора - 
размыкание ° и замыкание »/с помощью В опре- 
деляются как 

/оВ = (/ОВ)®В, 

/•В = (/®В)ОВ. 

Для иллюстрации геометрического поведения 
этих операторов полезно рассмотреть такие дву- 
мерные множества, как множество и структури- 
рующий элемент В, показанные в верхней части 
рис. 1. Этот рисунок иллюстрирует, что центр 
структурирующего элемента располагается на гра- 
нице преобразуемого множества; эрозия приводит 
к уменьшению множества X, а наращение - к его 
увеличению. Размыкание подавляет острые высту- 
пы и прорезает узкие перешейки в X, тогда как за- 
мыкание заполнят узкие заливы и малые отвер- 
стия, и таким образом ХОВ<^Хс:Х»В. Следователь- 
но, если структурирующий элемент В имеет регу- 
лярную форму, размыкание и замыкание можно 
рассматривать как нелинейную фильтрацию, кото- 
рая сглаживает контуры входного сигнала. Ясно, 
что форма и размер структурирующего элемента 
определяют природу и степень сглаживания. 

Таким образом, алгоритм реализации морфоло- 
гических операций для полутонового изображения 
заключается в следующем. Создается массив, тако- 
го же типа как изображение, который имеет такие 
же размеры как исходное изображение. В него за- 
писывается результат морфологической операции. 
Выбирается окно размером структурирующего эл- 
емента. Окно перемещается попиксельно по всему 
изображению и для точки оказавшейся в центре 
окна выбирается минимальное или максимальное 
(в зависимости от выполняемой операции, эрозия 
или наращение) значение из значений пикселей 
находящихся в пределах окна и это значение запи- 
сывается в массив. 
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В реальных технических системах, в которых 
решаются такие задачи как, идентификация но- 
мерных знаков автомобилей, распознавание тек- 
ста, анализ текстуры ткани и топологии печатных 
плат, возникает необходимость в обработке бинар- 
ных изображений. В подобных системах требуется 
высокая точность и скорость реализации алгорит- 
мов при минимальных затратах памяти. Тради- 
ционные методы морфологических преобразова- 
ний в данном случае являются не эффективными. 

Так как бинарные изображения представлены 
только двумя градациями яркости, например, 0 и 1, 
то выбор минимального или максимального эл- 
емента, в пределах окна структурирующего эл- 
емента, для каждой точки изображения является 
неэффективным. В этой связи нами предлагается 
модификация алгоритма выполнения базовых 
морфологических операций, эрозии и наращения. 
Полагаем, что форме принадлежат точки с интен- 
сивностью 1 . Структурирующий элемент выбираем 
размером ЬхЬ (структурирующий элемент может 
быть выбран любой формы, в данном случае для 
упрощения рассмотрения он выбран в форме ква- 
драта), хи у — координаты точки сканирования, ТѴх 
- ширина изображения, і(х,у) - значение интен- 
сивности в точке (х,у), ТѴу - высота изображения. 
Алгоритм имеет следующие шаги: 

1 . создается массив, такого же типа как изображе- 
ние, размером ЫххЫу, в котором будет содер- 
жаться результат операции; 

2 . координатам хи у присваиваем по значению Ь/2 ; 

3. вычисляем значение і(х,у); 

4. если і(х,у)= 1, тогда вычисляем по очереди зна- 
чение интенсивностей следующих точек: 
(х-1 ,у), (х— 1 ,у— 1 ) , (х,у- 1), (х+1,у-1), (х+1,у), 
(х+1,у+1), (х,у+1), (х-1,у+1); 

5. если і(х, у) =0, то идем на шаг 7; 

6. если одно из значений интенсивностей, вычи- 
сленных на шаге 4, равняется 0, то центр окна 
ЬхЬ располагается в точке (х,у) и всем точкам в 
массиве, попадающим в это окно, если это опе- 
рация эрозия присваивается значение 0, или, 
если это операция наращения 1; 

7. точке ( х,у ) в массиве присваиваем значение 0; 

8. если х<Ш-Ъ/ 2 и у<Иу-Ъ/ 2, тогда выполняем 
х=х+\ и идем на шаг 2; 

9. если х=Ш-Ъ/ 2 и у<Ыу-Ъ/ 2, тогда выполняем 
х=Ь/ 2, у=у+ 1 и идем на шаг 2; 

10. конец алгоритма. 

Оценка скорости реализации алгоритма пока- 
зывает, что количество операций присвоения точке 
значения интенсивности будет пропорционально 
лишь количеству точек в форме умноженное на 
размер структурирующего элемента. Отличитель- 
ной особенностью алгоритма является тот факт, 
что в нем исключены операции поиска минималь- 
ного или максимального элемента внутри структу- 
рирующего элемента. 


Далее приводится описание модификации ал- 
горитма реализации морфологических операций 
размыкания и замыкания для бинарных изображе- 
ний. Так, как смысл морфологически операций 
размыкания и замыкания заключается в том, что 
размыкание подавляет острые выступы и прорезает 
узкие перешейки в 2Г(рис. 1), а замыкание заполня- 
ет узкие заливы и малые отверстия, то предлагает- 
ся следующая модификация алгоритма выполне- 
ния операций размыкания и замыкания для бинар- 
ных изображений. Рассмотрим алгоритм реализа- 
ции операции замыкания по шагам (для операции 
размыкания выполняется то же самое, что и для за- 
мыкания, только вместо нуля единица, а вместо 
единицы ноль). 

Выбираем те же обозначения, что и для базовых 
морфологических операций, ^-счетчик. 

1 . Координатам хи у присваиваем по значению 0, 
^0; 

2. вычисляем значение і(х,у); 

3. если і(х,у)= 1, идем на шаг 9; 

4. если і(х,у)= 0, вычисляем /(х+1 ,у); 

5. если /(х+1,у)=0, выполняем к=к+ 1, х=х+1 и вы- 
числяем значение Нус, у) до тех пор, пока і(х,у) не 
станет равным 1 или пока х не станет равным /Ѵх; 

6. если х=/Ѵх, тогда ^=0 и идем на шаг 10; 

7. если І(х,у)= I и к<=Ь, тогда всем точкам от 
(х-к,у) до (х,у) присваиваем значение 1; 

8. к= 0; 

9. если х<Ах и у<Лф, х=х+ 1 и идем на шаг 2; 

10. если х=/Ѵх и у<Лф, х=0, у=у+ 1 и идем на шаг 2; 

11. координатам хи у присваиваем по значению 0, 
^0; 

12. вычисляем значение і(х,у); 

13. если і(х,у)= 0, идем на шаг 19; 

14. если і(х,у)= 1, вычисляем /(х,_у+1); 

15. если /(х,у+1)=0, выполняем к=к+ 1, у=у+ 1, и 
вычисляем і(х,у) до тех пор, пока і(х,у) не станет 
равным 1 или пока у не станет равным А !у; 

16. если у=Му, тогда к= 0 и идем на шаг 20; 

17. если І(х,у)= I и к<=Ь, тогда всем точкам от 
(х,у-к) до (х,у) присваиваем значение 1; 

18. к=0\ 

19. если у<Му и х<Ах, у=у+ 1 и идем на шаг 12; 

20. если у=Иу и х<7Ѵх, х=0, х=х+ 1 и идем на шаг 12; 

21. конец алгоритма. 

Оценка скорости реализации алгоритма пока- 
зывает, что количество операций присвоения точке 
значения интенсивности будет пропорционально 
лишь количеству точек одинаковой интенсивности 
располагающихся в промежутке, размером равным 
или меньшем, чем размер структурирующего эл- 
емента, между точками другой интенсивности. От- 
личительной особенностью алгоритма является тот 
факт, что в нем исключены операции поиска мини- 
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Технические науки 




Рис. 2. Применение быстрых алгоритмов операций замыкания и размыкания: а) исходный сигнал, б) замыкание структуриру- 
ющим элементом 5х 7, в, г) размыкание структурирующим элементом 1x3 и 15х 1 


мального или максимального элемента внутри 
структурирующего элемента. А также исключена 
повторная замена точек форм, которая бы осу- 
ществлялась при поочередном использовании мо- 
дифицированных операций эрозии и наращения. 

Описанные выше быстрые алгоритмы опера- 
ций замыкания и размыкания применяются для 
предварительной обработки изображений автомо- 
биля, в задаче обнаружения номерной пластины. 
На рис. 2 представлено последовательное примене- 
ние набора морфологических операций к бинари- 
зованному изображению. На рис. 2, б, представлен 
результат применения операции замыкание струк- 
турирующим элементом размером 5x1. Размер 
структурирующего элемента равняется расстоянию 
между буквами на номерном знаке автомобиля, т. е. 
размеры номерной пластины должны быть заранее 
известны. 

Результат для наглядности представлен черным 
и серым цветом, т. е. точки с интенсивностью 0 и 
150 принадлежат одной форме. 

На рис. 2, в представлен результат операции 
размыкания, примененной к изображению рис. 2, 
б, размер структурирующего элемента выбран раз- 


мером меньшим, чем высота букв 1x3. На рис. 2, г 
- результат операции размыкания, примененной к 
изображению рис. 2, в, размер структурирующего 
элемента выбран размером меньшим, чем длина 
номерной пластины 15x1. 

Данный метод, обладает следующими недостат- 
ками. Это проблема выбора порога бинаризации. 
Данную проблему, возможно, решить с помощью 
алгоритмов адаптивной бинаризации. Второй не- 
достаток связан с трудностью определения разме- 
ров структурирующего элемента, т. к. он связан с 
размерами символов. 

Программа была разработана и реализована в 
визуальной среде программирования Ое1рІіі7. Во 
время проведения испытаний, на компьютере с 
конфигурацией А1ѴГО АіЫоп 2, 1 ГГц и оперативной 
памятью объемом 256 Мб, для изображения разме- 
ром 1024x768 время исполнения составило 20 мс. 
Время реализации морфологических операций 
традиционным путем составило 900 мс. 

Отличительной особенностью предложенных 
алгоритмов является исключение операции поиска 
минимального или максимального элемента вну- 
три структурирующего элемента, что существенно 
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уменьшает время реализации данных алгоритмов 
для бинарных изображений. Для операций замыка- 
ния и размыкания исключена повторная замена 
точек форм, которая бы осуществлялась при по- 
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